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Abstract 

The dilute A3 model is a solvable IRF (interaction round a face) model with 
three local states and adjacency conditions encoded by the Dynkin diagram 
of the Lie algebra A3. It can be regarded as a solvable version of an Ising 
model at the critical temperature in a magnetic field. One therefore expects 
the scaling limit to be governed by Zamolodchikov's integrable perturbation 
of the c = 1/2 conformal field theory. Indeed, a recent thermodynamic Bethe 
Ansatz approach succeeded to unveil the corresponding Eg structure under 
certain assumptions on the nature of the Bethe Ansatz solutions. In order to 
check these conjectures, we perform a detailed numerical investigation of the 
solutions of the Bethe Ansatz equations for the critical and off-critical model. 
Scaling functions for the ground-state corrections and for the lowest spectral 
gaps are obtained, which give very precise numerical results for the lowest 
mass ratios in the massive scaling limit. While these agree perfectly with the 
Eg mass ratios, we observe one state which seems to violate the assumptions 
underlying the thermodynamic Bethe Ansatz calculation. We also analyze the 
critical spectrum of the dilute A3 model, which exhibits massive excitations 
on top of the massless states of the Ising conformal field theory. 

05.50.+q, ll.25.Hf, 75.10.Hk 
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I. INTRODUCTION 



The Ising model |IJ is, without doubt, one of the most frequently studied and best 
understood lattice models of classical statistical mechanics. Although Onsager's solution 
of the two-dimensional Ising model without external field dates back half a century already, 
no analytic solution of the 2D Ising model in a magnetic field has been found. 

However, the situation is somewhat different if considered from the viewpoint of field 
theory. The critical 2D Ising model corresponds to the c = 1/2 conformal field theory 
(CFT) of a massless Majorana fermion. It was Zamolodchikov |7]-|9| who noticed that 
a (symmetry-breaking) perturbation of this CFT with the relevant spin density operator 
(which has conformal dimensions (A, A) = (1/16,1/16)) preserves infinitely many conser- 
vation laws and therefore leads to an integrable quantum field theory. The corresponding 
minimal theory contains eight massive particles with factorized (purely elastic) scatter- 
ing [|TX?|^r5j] ; the particle masses and the S'-matrix elements are related to the exceptional Lie 
algebra E 8 . 

This integrable field theory describes an appropriate scaling limit of the 2D Ising model 
in a magnetic field. Numerically, the predictions for the lowest mass ratios have been 
verified by several authors ||16|-|20||. As these results rely on relatively small size transfer 
matrix calculations or on Monte Carlo simulations, they provide rather crude checks for 
the lowest mass ratios only. Furthermore, the larger masses cannot be obtained directly as 
these lie above the two-particle threshold of the lightest particles and hence are buried in 
the continuum of scattering states. However, one can extract indirect information about 
the higher masses and the scattering amplitudes from the finite-size behavior of the lower 
Another, rather different, approach employs a truncated conformal Hilbert 
which however does not make the connection to the lattice model. 
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Similar theoretical and numerical investigations have also been performed for a variety of 
lattice models such as the quantum Ising chain |16| , the Lee- Yang model [0] , the three-state 
Potts model [T24|], the tricritical Ising respectively Blume-Capel model [5S|-27], the Ashkin- 
the Z(N) models p9| , and the integrable RSOS models |j30| . Recently, 



Teller model 28 



this has been extended to the study of form factors and correlation functions |3~T|-|3l)|| and 
non-integrable perturbations |3(| . 

The discovery of a new class of lattice models pT|-piS[| which are solvable in the presence 
of a symmetry-breaking field [38 has changed the situation considerably. Although we still 
cannot solve the 2D Ising model in a magnetic field, we now know a solvable model (the 
so-called dilute A 3 model) which belongs to the same universality class: it has a critical 
point of Ising type [p8| , ^0| , and it is solvable in a symmetry-breaking field that corresponds 
to the spin density perturbation of the CFT Therefore, one expects this model to show 
the same properties in the scaling limit as the non-integrable Ising model in a magnetic field. 
Note, however, that the two models are rather different lattice models: the dilute A3 model 
contains three local states, whereas the Ising spins take two values only. 

The dilute A^ models are solvable by Bethe Ansatz (BA) By a thermodynamic 

Bethe Ansatz (TBA) approach ||42|| , the mass ratios and S'-matrix of the E 8 field theory 
have been obtained. However, this approach relies on assumptions on the nature of the 
BA solutions, which could not be substantiated by numerical solution of the BA equations 
(BAE) at criticality. Connections between the A 3 model and the exceptional Lie algebra E 8 
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are also revealed by the existence of corresponding Roger s-Ramanuj an identities [44 . 

For our numerical investigation of the off-critical BAE, our main motivation was to 
check the assumptions on the conjectured string structure of the BA solutions. This could 
be achieved for the states with lowest masses, yielding at the same time very accurate 
numerical results not alone for the mass ratios, but also for the complete scaling functions. 
However, we also found one discrepancy with the predictions of Ref. ||42|| . Within the range 
of scaling parameter values considered, the characteristic string type for a particle of mass 
m 4 does not match the proposed structure. 

As a by-product of our investigation, we found that the critical spectrum also contains 
massive excitations besides the well-known conformal spectrum. This is one property that 
distinguishes the dilute A3 model from the Ising model which does not have such excitations. 
However, this does not contradict universality as that only concerns the universal properties 
of the systems in the vicinity of the critical point, which are the same for both models. 
The BA solutions corresponding to the massive excitations show the same type of strings as 
those conjectured for the massive scaling limit. Already for relatively small systems, these 
solutions lie very close to singularities of the BAE, which makes the numerical treatment 
difficult. Estimates of the corresponding mass ratios, which show no resemblance of the E§ 
structure, are presented. 

This paper is organized as follows. In Section [TT|, the dilute A 3 models is introduced and 
the BAE are presented. The critical spectrum predicted by CFT and the behavior in the 
massive scaling limit are described. Furthermore, the assumptions on the form of the BA 



solutions formulated in Ref. [421 are summarized. 



Section |TTTj contains the results of our numerical investigations. After a brief description 
of our numerical approach, we commence with the critical case, where not only the scaling 
dimensions predicted by CFT are considered, but also additional massive excitations in the 
spectrum are identified. Then, we show how the BA solutions behave as one goes away from 
criticality. Finally, we present our numerical scaling functions for the lowest gaps and the 
numerical results for the mass ratios in the massive scaling limit. 



II. THE DILUTE A, MODEL 



Because an immediate specialization to the A3 case does not considerably simplify the 
equations, we keep the first part of this section on a more general level and consider the 
dilute A^ models. This also serves us to avoid an unnecessary repetition in a forthcoming 
publication f4"5fl . 

The dilute A^ models are IRF (interaction round a face) models HH] of RSOS (restricted 



solid-on-solid) type |47j on the square lattice with adjacency conditions encoded by the 
Dynkin diagram of the Lie algebra A^. In contrast to the usual (non-dilute) RSOS models 
built on this adjacency graph (these are the ABF models 10), the effective adjacency graph 
of the dilute A L models contains loops which connect each node to itself, see Fig. [j]. In 
other words, one considers configurations of L local states or "heights" (labeled 1,2, ... ,L) 
on the vertices of the square lattice subject to the additional requirement that states on 
neighboring (adjacent) lattice sites may differ at most by one. One then defines a statistical 
model by assigning Boltzmann weights to the elementary plaquettes (faces) of the lattice. 
The face weights of the dilute A^ models are 
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Here, a = 1, 2, ... L labels the heights, and the possible values of the variable A are deter- 
mined by L 



. The crossing factors S a are given by 

, ^i(4aA) 
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and i?4(u) are standard theta functions of nome q with |g| < 1 f48| . These face 



weights satisfy the Yang-Baxter equation and therefore lead to an integrable lattice model. 



The corresponding row transfer matrices |5| for fixed value of q form a one-parameter 
commuting family in the spectral parameter u. 

In general, the dilute models have four different branches for any value of L 
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distinguished by two possible values of A with two regimes each for the spectral parameter 
u. Here, we are only interested in the dilute A3 model (L = 3) in a particular branch 
P^ , [4D|J42]| where the model behaves like the Ising model in a magnetic field. For this case, 
A = 57r/16 and the spectral parameter u lies in the interval < u < 3A. It is the parameter 
q that acts like the magnetic field, and the Ising critical point corresponds to q = 0. 
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A. Bethe Ansatz Equations 



The eigenvalues A(u) of the row transfer matrix for a system of size N with periodic 
boundary conditions have the following form P2| , pE3 | 



0i (u - 2A) x (u ~ 3\)\ N A(u - Uj + A) 



1 0i(2A)0i(3A) J LL&^u-Uj-X) 

I 0i(u) 0i(u — 3A) \ N rr 0i (w — Uj) 0i(m — Uj — 3 A) 



0i (2A) 0i (3A) J Al ^ ( u _ Uj . _ A) 0! (u - Uj - 2A) 

i / 0i(u)0i(u-A) \ JV Jt ty 1 ( M - Mj --4A) 
V 0i(2A)0i(3A) J ^0i(u-Uj-2A) 1 ' J 

where w = exp(i7r/ (L + 1)), and where the Uj, j = 1, . . . , N, form a solution of the set of iV 
coupled BAE 



0i {uj - A) A N _ yr M% ~ M fc ~ 2A ) - ^fc + A) , . 

0i(uj + A)j ~ ll0i(uj-u fe + 2A)0i(uj-u ifc -A) 1 ' J 



Here, £ = 1, . . . , L labels a sector related to the braid limit eigenvalues of the row transfer 
matrix at criticality (q = 0). Whereas it is believed to be true that all eigenvalues of the 
transfer matrix are of the form ( |2.3| ), the converse is certainly wrong (at least off-criticality) 
- the BAE ( |2.4| ) in general allow for many additional solutions which do not correspond to 
proper eigenvalues of the transfer matrix. For any numerically found solution of Eq. ( |2.4j ) one 
thus has to check whether it gives a proper eigenvalue of the transfer matrix (see Section |nj 
below). 



The theta function 0i(w) shows the (quasi-)periodicity properties |48 



1 (u + 7 r) = -0i(u) = 0i(-u) (2.5a) 

■& x (u + wr) = --e-^dAu) (2.5b) 
Q 

where q = exp^inr) (0 < q < 1) with r G iR. Consider a solution {ui, . . . , u^} of the BAE 
(|2.4|) in sector £. Clearly, nothing is changed if a multiple of 7r is added to any of the roots 

Uj. 

More interesting is the result of adding rnir (n G Z) to a root Uj. In this case, one obtains 
a solution of (|2.4| ) with u e replaced by exp(4mA)a; . Moreover, from Eq. ( |2.3|) it is evident 
that these two solutions correspond to the same eigenvalue. This means that the "sectors" 
£ — 1, . . . , L lose their significance in the off-critical case (q ^ 0) since we can always adjust 
£ by adding or subtracting suitable multiples of ttt to some of the roots. Of course, this also 
yields many possibilities in which the phase factors cancel, and hence Eq. ( [2.4D is recovered 
without any alteration, which means that a single solution of the BAE can be presented in 
many ways differing by addition and subtraction of suitable multiples of tit to some of the 
roots. 

For the case of interest, it turns out that the largest eigenvalue A (u) of the transfer ma- 
trix is given by a purely imaginary solution of the BAE (|2.4j) in the sector £ = 1. Therefore, 
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in place of the roots Uj, we prefer to use Vj = iuj in what follows; thus the largest eigenvalue 
corresponds to a set of real BA roots Vj. General solutions to the BAE will however involve 
complex roots, which in the large N limit typically arrange into so-called strings (subsets of 
roots with approximately the same real part, in terms of the Vj). This is shown schematically 
in Fig. 0, where the horizontal line represents the real axis. Note that the "string content" 
of a particular solution might well depend on the variable q, an example of this behavior 
will be given in Section |111 C| below. 



B. Conformal Spectrum at Criticality 

The dilute A 3 model has a critical point of Ising type at q — 0. Thus its critical limit 
(given by q = 0, iV — ► oo) is described by the c = 1/2 CFT with scaling dimensions 
A e {0, 1/16, 1/2}. This has drastic consequences for the spectrum of the transfer matrix 
in the critical limit. Consider the scaled spectral gaps 

*j = ^ log(A( 0) /Af ) , (2.6) 

where are the eigenvalues of the transfer matrix with periodic boundary conditions at 

the isotropic point u = 3 A/2, A ^ being the largest eigenvalue. In the critical limit, the 
spectrum of scaling dimensions Xj consists of "conformal towers" of states labeled by pairs 
( A + r, A + r) with scaled energy Xj = A + r + A + r, conformal spin Sj = A — A and 
momentum pj = r — r, where r,r G No- They form representations of two commuting 
Virasoro algebras with central charge c = 1/2. The degeneracies of the descendant states 
can be read off from the character functions of their irreducible representations with highest 
weight A. These can be written in the form j|9 



x (z) = ;>> 4 " 2+ "n y (^) 

= l + z 2 + z 3 + 2z A + 2z 5 + ... (2.7a) 



X±{z 
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= zi (1 + z + z 2 + z 6 + 2z* + ...) (2.7b) 

oo 

= z^l[(l + z m ) 

16 ■ L 

m=l 

= zTe (1 + z + z 2 + 2z 3 + 2z 4 + . . .) (2.7c) 
where IL v (z) is the generating function of the number of partitions 

oo 

n v^) = II T3^r • ( 2 - 8 ) 

m=l 

Moreover, the central charge c = 1/2 manifests itself in the finite-size corrections of the 



largest eigenvalue [50 
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where /o denotes the bulk free energy for u = 3 A/2 
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(2.9) 



C. Scaling Limit 

Taking into consideration the parameter q, one can approach the critical point by simul- 
taneously performing the two limits g — f and N — > oo keeping the scaling variable 



Q 



N 15/8 



(2.10) 



constant. Here, /i = corresponds to the critical limit discussed above, fi — > oo to the 
massive limit where Zamolodchikov's E 8 field theory results apply. In the scaling limit, the 
appropriately scaled spectral gaps 



F 3 = MAo/A; 



(2.11) 



become functions of the scaling variable \i alone. For the largest eigenvalue, we can also 
define a scaling function 



q' 8 / 15 [log(A ) + N ■ f 



(2.12) 



for the finite-size corrections of the largest eigenvalue, where fo = liniAr^oof— log(A )/A] is 
the bulk free energy again. Note that the face weights ( |2.1| ) for nome q and — q are related 
by symmetry, thus it suffices to consider positive values of the nome. 
In the massive limit (fi — > oo), the ratios 



Rj 



F 1 



(2.13) 



approach the particle mass ratios of the corresponding massive field theory. The masses of 
the eight stable particles are proportional to the entries of the Perron-Frobenius eigenvector 
of the Cartan matrix of the Lie algebra E 8 and their ratios (ordered by magnitude) are given 
by 



mi/m 


= 1 






m2 /m 


= 2 


cos(7r/5) 


= 1.618 034... 
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cos(7r/30) 
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m^/m 


= 4 


cos(7r/5) cos(77r/30) 


= 2.404 867... 


m 5 /m 


= 4 


cos(7r/5) cos(27r/15) 


= 2.956 295... 


m e /m 


= 4 


cos(7r/5) cos(7t/30) 


= 3.218 340... 


m-jjm 


= 8 


cos 2 (tt/5) cos(7tt/30) 


= 3.891 157... 


mg/m 


= 8 


cos 2 (7r/5) cos(27t/15) 


= 4.783 386... 



(2.14) 



where m = mi defines the mass scale. The Dynkin diagram of Eg is shown in Fig. || where 
our labeling of nodes follows the ordering of the masses. 
In the TBA calculation of Ref. 



42], each of the corresponding eight massive particles is 



associated to a particular string in the BA solution. In addition, it is assumed that, apart 
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from these eight strings and the one-strings forming the vacuum solutions, no other string 
types occur in thermodynamically relevant quantities. With these assumptions, the TBA 
equations imply that the hole-type excitations in the one-strings have vanishing density. 
The density of one-strings can be eliminated from the calculation, which finally leads to the 
8x8 scattering matrix of the Eg factorized scattering theory p2 |. 



In Table |, the nine string types are given, labeled by t — 0, 1, ... 8 in the order used 
in Ref. | j42|| , which corresponds to the usual labeling of the Eg Dynkin diagram rather than 
that used in Fig. |^. From the data of Table p], the string of type t consists of roots of 
the form 



+ |(A<" + 164") (2.15) 
where denotes the center of the string on the real line. 



„(«) = „<*> , 
k 32 



III. NUMERICAL RESULTS 



For the numerical treatment of the BAE ( |2.4|) we used a modified Newton method |5T 
The calculations were performed in extended precision FORTRAN (with 16 byte real num- 
bers) on an IBM workstation. 

As a first step, we solved the BAE at criticality (q = 0) for small system size and 
compared the corresponding eigenvalues ( |2.3j ) to those obtained by direct diagonalization 
of the transfer matrix using high precision arithmetics. This also served as a test of the 
performance of our programs. In particular, the complete set of solutions corresponding to 
zero-momentum eigenstates of the transfer matrix for systems of size N < 6 was obtained. In 
this way, the structure of the solutions for the largest eigenvalues could be identified. These 
solutions were followed as a function of the elliptic nome q. The results of this analysis were 
then generalized to larger system size by looking for solutions of the same type which differ 
by additional real roots only. 

Since we are mainly interested in the mass ratios in the scaling limit, we exclusively 
considered zero momentum states. For solutions of the BAE, the momentum is obtained 
from the eigenvalue A(0) (|2.3|) at spectral parameter u = 0, where the transfer matrix 
reduces to a shift operator. 

Because we want to follow solutions for varying q, which means that we have to repeat 
the calculation for a large number of values of q, we limited ourselves to systems of size 
N < 100. However, in certain cases the numerical accuracy of the calculations becomes 
the major problem. This happens for instance for BA solutions which contain roots that 
are very close to singularities of the BAE, or for solutions where differences between roots 
(respectively their real parts) become extremely small. These situations typically show up if 
one considers large values of the elliptic nome q (where the meaning of "large" depends on the 
system size), but also for the massive excitations of the critical dilute A 3 model discussed in 
Section [U1B| below. As far as we can see, there are only two ways to go beyond the limitations 
imposed by these numerical problems; either by using higher precision arithmetics, or, to a 
lesser extent, by adjusting the program according to each specific type of BA solution (which 
we did to some extent by treating the ubiquitous complex conjugate pairs of BA roots with 
imaginary parts close to ±ll7r/32 in a special way). 
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A. Critical Conformal Spectrum 



In the critical limit (q = 0, N — > oo), the largest eigenvalues of the transfer matrix 
organize according to the characters of representations of the Virasoro algebra with central 
charge c = 1/2 as described in Section |II B| above. Schematically, the resulting values for 
the scaling dimensions x in the zero-momentum sector and their distribution into the three 
sectors labeled by £ = 1, 2, 3 are shown in Fig. f|. 

In Table [TJ, numerical values for the central charge and the lowest scaling dimensions 
obtained from the corresponding solutions of the BAE for systems of different sizes are 
presented. Here, Eq. ( |2.9| ) with the exact value of the bulk free energy |52]J5B[ ] of the dilute 
A3 model was used. All zero-momentum states with a critical scaling dimension of x < 7 have 
been found, compare Fig. f|, except one of the four excitations with conformal dimensions 
(1/16 + 3,1/16 + 3) for which, despite some effort, we have not been able to find the 
corresponding solution of the BAE. In particular, this includes all cases without degeneracy, 
i.e., where only a single zero- momentum state of a certain scaling dimension occurs in the 
conformal tower. 

To keep the three remaining BA solutions with scaling dimension xg = 6 + 1/8 apart, 
we distinguish them by subscripts a, b, and c. Actually, it turns out that the two solutions 
denoted by b and c yield the same eigenvalue even for finite systems, which remains true also 
off-criticality. Thus, we do not need to consider them separately unless we are interested in 
the BA solutions, which of course are different, see Section [111 C| below. 

While the numerical data are in perfect agreement with the exact values, it should be 
noted that the string structure of the critical BA solutions have no apparent similarity to 
those proposed in Ref. E3. We shall come back to this point later. 



B. Massive Excitations at Criticality 

Analyzing the spectrum of largest eigenvalues of the transfer matrix, one observes that, 
in addition to the conformal spectrum discussed above, the spectrum contains massive ex- 
citations. More precisely, these are states with eigenvalues for which the quantities 

Vk = log(Ai 0) /Ai 0) ) (3.1) 

converge to a non-zero limit (mass) as A" — > 00, compare Eq. fl2.6|) for the conformal states 
where this limit gives zero. 

Apparently, a number of different masses is involved, which are again characterized 
by particular strings in the BA solutions. Strikingly, the strings which we observed are 
among those listed in Table |. We have found explicit solutions containing strings of type 
t = 1, 2, 3, 4, and 7. In all cases, the roots forming these strings are — even for small systems 
- located extremely close to the singularity of the BAE which they approach for A" — > 00. 
This asks for a careful numerical treatment, and severely limits the system sizes we can 
treat. With the numerical accuracy of our program, we were able to get reliable results for 
systems of at most 20-30 sites, depending on the particular solution. 

An example is given in Table [ITT]. Here, the solutions obtained by adding a number 
of "massive" two-strings (t = 1) to the ground-state solution of a system of size A^ = 4 
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are presented. From the eigenvalues for larger systems, one can see that each additional 
"massive" two-string just adds the same mass, thus substantiating our interpretation in 
terms of massive particles. 

Apparently, it is possible to add arbitrary numbers of these "massive" strings not only to 
the ground-state, but to each solution corresponding to a conformal state. This means that 
the spectrum (at the critical point q — 0) contains infinitely many copies of the complete 
conformal spectrum, shifted with respect to the ground-state by a mass determined by the 
collection of "massive" strings. 

As to the observed "massive" strings, they fall into two classes: for t = 1,3, and 7, we have 
a single solution of the BAE, whereas for t = 2 and 4 the solutions belong to eigenvalues 
which are doubly degenerate (for finite systems). Labeling the masses associated to the 
"massive" string type t by M t , we obtain the numerical estimates 

M 7 /M 1 « 1.85186 M 2 /M 1 « 1.78 , , 

M 3 /M 1 « 3.17213 M 4 /Mx « 3.78 ^ ' 

for the ratios of the masses with respect to the smallest mass which is associated to the 
two-string (t = 1). Note that these numbers are calculated from the eigenvalues of the 
transfer matrix with spectral parameter u = 3 A/2. 

This situation is similar to what is observed in the Hubbard model, where one has 
two types of excitations: so-called "holons" (carrying charge but no spin) and "spinons" 
(carrying spin but no charge) Generically, both excitations are massless, and the 

corresponding field theory consists of two coupled c = 1 CFT. However, if the filling fraction 
(number of electrons per lattice site) is chosen to be precisely one, the holons become massive 
while the spinons stay massless. At low energies, this theory is then described by a c = 1 



CFT with additional massive excitations [56] 



C. Bethe Ansatz Roots Off-Criticality 

As mentioned above, the string structure of the BA roots at criticality does not seem 



to agree with the predictions of Ref. ||42|| . However, if one follows a particular solution as 
a function of the nome q, one finds that in many cases one encounters singularities of the 
BAE where the string structure changes. 

The simplest example is given by the second largest eigenvalue, which corresponds to 
the conformal spin density field with dimensions (1/16,1/16) in the critical limit and to 
the lightest massive particle in the massive scaling limit. At criticality, the corresponding 
solution of the BAE differs from that of the largest eigenvalue (which contains real roots v j 
only) by a single root with imaginary part tt/2. However, as can been seen from Fig. [|, as 
q is increased from q = 0, the complex root and one real root approach each other, until 
their real parts agree. At that point, the real parts of the two roots stay the same, but they 
move in the imaginary direction until they form a complex conjugate pair of roots (hence a 
two-string), which then persists for large values of q. 

Comparing systems of different size, this mechanism stays always the same, the additional 
real roots just play the role of spectators. Moreover, the transition from one string type 
to another takes place at approximately the same value of the scaling parameter \i (|2.10| ) 
as the system size is increased. More precisely, the scaling parameter values approach a 
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non-zero limit as the system size tends to infinity. This implies that although q — > in the 
scaling limit, for large values of [i (and in particular for the massive scaling limit \i — > oo), 
the two-string is the relevant type of solution. 

The same scenario applies to the other excitations we considered. In Table |IV|, we 
compiled numerical values of the non-real BA roots for systems of size N = 100, both at 
criticality (q = 0, i.e. \l — 0) and for a rather large value /i ~ 39.4 of the scaling parameter 
(corresponding to q — 7/1000 for iV = 100). Note that for q ^ the real part of the roots Vj 
can be shifted by multiples of — mr = — log(g) G JR+ (0 < q < 1) as discussed in Section [II A 
above. Using this property, we arranged the solutions such that the values of the non-real 
roots have a real part in the interval [0, — log(g)). Actually, as Table [TV] suggests, all non- 
real roots are located in the vicinity of — log(g)/2 for sufficiently large scaling parameters \i. 
This also distinguishes these solutions from the "massive" strings at criticality considered 
above, whose centers are located close to the origin. Furthermore, the values of the nome 
q where changes in the string patterns for N = 100 occur are also indicated in Table [IV]. 
All five single-particle states show at least one such point, whereas some of the two-particle 
states do not change as q is varied. 

Although the relation between the scaling exponents on one end and the content of 
massive particles on the other is determined by the field theory, we are not aware that it 
has been calculated. Our results given in Table |VJ give an unambiguous connection for the 
lowest states, but they do not suggest an obvious pattern that can be generalized to the 
higher excitations. 

Comparing with the proposed string types of Ref. given in Table [I] and by Eq. ( p,15|) , 
one finds that the solutions for large scaling parameter agree with the predictions apart 
from the two states which contain a particle of mass m 4 . Here, in place of a five-string with 
imaginary parts ±7r/4, ±37r/8 and tt/2, the solution found numerically on first view looks 
more like a seven-string with imaginary parts ±7r/32, ±77r/32, ±137r/32 and tt/2. However, 
in comparison to the other single-particle states, the real parts of the roots forming this 
"string" are quite far apart, which might be an indication for another change of pattern 
occurring at a larger value of q, which of course we cannot exclude on the basis of our 
numerical data. 

Let us have a closer look at the single particle state with mass m 4 . Figures |] and [7] 
show the BA roots of the corresponding solutions for systems with N = 10 and = 20 
sites, respectively. For convenience, we plotted the real and imaginary parts of the roots, 
normalized by — log(g) and by tt as in Fig. [5], against /j, 8 ^ 15 = q 8 / 15 N. There are two obvious 
points where non-analyticities can be seen in the imaginary parts, which both are situated 
at approximately the same value of the scaling parameter for the two different systems. 
However, Fig. |6] may suggest that another change will occur at a larger value of fi than 
numerically accessible by our routines, as the real parts of one real root and one of the two- 
strings approach each other. But, comparing with the same region of Fig. [7], no indication 
of this behavior remains, and the same holds true for the larger systems (up to N = 100) 
which we examined. From this, we conclude that either there is no further change in the 
string pattern, or it has to occur for a very large value of the scaling parameter /z, which 
however seems rather unlikely to us. 

Though this result is somewhat inconclusive, let us examine the string pattern of this 
solution in more detail. In particular, we are interested in the question whether there are 
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holes in the one-strings in this case. To see this, we consider the "phase function" 



<p{v) = ± lOf 



_ -e f M iv + A) ' 



N 



N 



T-r ~&i(iv — ivk — 2A) — ivk + A) 
^ - iv k + 2A) - iv k - A) 



(3.3) 



which is basically the logarithm of the BAE ( |2.4| ), but now v k denote the roots of our 
particular BA solution, and we consider tp as a function of the complex variable v. By 
definition, <p(vj) G Z for all the roots Vj, j — 1, . . . , N. Restricting to real values of v, tp{v) 
therefore takes integer values at all real solution v = Vj G R, hence for all one-strings. 

Fig. § shows the <p(v) for the solution under consideration. Here, the size of the system 
is N = 20, and q = 1/5. The horizontal lines are drawn at integer values (the precise 
numbers are not relevant as they depend on the choice of branch in Eq. ( |3.3| )), and the 
crosses located on intersections of these lines with the graph of ip(v) denote the position of 
the real roots. The vertical lines are placed at the real parts of the remaining seven non-real 
roots. As observed above, the single root with imaginary part tt/2 lies somewhat separated 
from the three complex conjugate pairs which are very close together (such that the three 
vertical lines appear as one thicker line in the figure). Moreover, the location of this single 
root coincides precisely with a hole in the one-strings, which may seem to disagree with the 
findings of Ref. ||42|| . However, the behavior of the roots as a function of parameters like 
q suggests that the roots with imaginary part tt/2 are simply an alternative locus of the 
one-string. 

On the basis of the numerical data, it thus seems that the string structure associated 
to the mass m 4 consists of the rather complicated pattern of a seven-string involving three 
complex conjugate pairs of roots (with imaginary parts ±7r/32, ±77r/32 and ±137r/32) and 
a neighboring root at in/ 2 which comes together with a hole in the one-string distribution. 
Of course, in the infinite size limit the distance between the real parts of these roots becomes 
infinitesimal. It is not obvious what happens in the TBA calculation of Ref. j|2) if their 



string type t — 8 (see Table |) is replaced by the observed pattern, especially as we did 
not see the strings associated to the three largest masses. Nevertheless, we believe that the 
main ideas underlying the TBA treatment in Ref. [H2J are correct, and that only details of 
the calculation would be affected. 

Before we move on to the discussion of the scaling function, two short remarks regarding 
the BA solutions are in order. It should be noted that the "massive" strings mentioned 



previously in Section [D [B|, which yield the massive excitations at criticality, do not undergo 
similar changes. These strings are always located very close to singularities of the BAE, and 
move even closer as the nome q is increased. Moreover, their centers lie at the origin, see 
e.g. Table |TJ, whereas for the states belonging to the E 8 integrable field theory the strings 
cluster around the value — log(g)/2. Finally, let us mention that a similar behavior of BA 
roots has recently been observed in the solution of the XXZ Heisenberg chain, where the 



boundary twist acts as the varying parameter [57 
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D. Scaling Functions and Mass Ratios 



Let us now turn to the results for the scaling functions F (n) ( |2.12| ) and Fj(fi) fl2.11| ). For 
the latter, all zero-momentum states with a critical scaling dimension of x < 7 are considered, 
except the one missing excitation with conformal dimensions (1/16 + 3, 1/16 + 3). In all 
cases, we use data from systems of size N = 50, N = 75 and N = 100. For these sizes and 
the values of q we considered, it turns out that the corrections to scaling are so small that 
one can hardly recognize them in our figures. 

In Fig. H the (natural) logarithm of F (fj,) (which as defined in Eq. ( [2.12|) is positive) is 
shown as a function of /i 8//15 = q 8 ^ 15 N ( |2.10| ). Apart from the behavior for small q, the plot 
is nicely linear, showing the exponential decrease of Fq as a function of /i 8//15 down to about 
exp(— 60) where the difference between the eigenvalues for the finite sizes and the bulk limit 
value becomes smaller than our numerical precision, which also proves the performance of 
our numerical routines. Small deviations from scaling can be seen for larger values of q. 

The scaling functions for the excitations Fj(n) ( |2.11| ) are displayed in Fig. [10]. The curves 
shown are piecewise linear plots connecting data points obtained from systems with N = 100. 
For convenience, we again used /i 8/15 = g 8 / 15 iV (f2~T0| ) on the horizontal axis. Qualitatively, 
the scaling functions agree with the results of previous numerical calculations for the Ising 
model in a magnetic field fi7|-|r5f and with the results of the truncated fermionic space 
approach, compare the figures given in Ref. f23|. The one-particle states can be recognized 



by their characteristic minima |T7|] . This becomes clearer when we consider their ratios 
Rj(fj>) ( p.!3|) which are presented in Fig. [11. Here, we again show individual data points 
which were obtained from systems of size N = 50, N = 75, and iV = 100, which obviously 
were large enough to keep corrections to scaling small. The ratios Rj(fi), given here against 
the scaling variable fi, are labeled by the corresponding conformal dimensions at criticality. 
We also indicated the single-particle mass ratios ( |2.14j ) of the E 8 field theory. Clearly, the 
agreement is convincing, more detail on the approach of the massive scaling limit is contained 
in Table |V[ 



IV. CONCLUSIONS 

The spectrum of the dilute A 3 model has been studied by numerical solution of the 
Bethe Ansatz equations, both at and off-criticality. This gives the correspondence between 
the lowest states of the conformal spectrum of the critical Ising model and those of the 
massive E 8 field theory (note that the connection given in Ref. |T7| contains an obvious 
mistake). 

At criticality, the spectrum consists for one part of massless excitations described by 
c = 1/2 minimal CFT as in the critical Ising model, but in addition to these it contains 
massive excitations. Obviously, the masses of these excitations are linked to the appearance 
of particular strings in the Bethe Ansatz solutions, which can be interpreted as massive 
particles. In our numerical analysis, we have seen at least five such "massive" strings. It 
is intriguing that all these strings (showing up in solutions at criticality) are among the 



conjectured list of "thermodynamically significant string types" in the TBA analysis |42 
of the massive scaling limit. This is strongly suggestive of a connection between the Lie 
algebra E 8 and the masses of these excitations. 
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However, these are not the states which yield Zamolodchikov's Eg field theory of the Ising 
model in a magnetic field. Those correspond to massless excitations at the critical point, 
which develop a mass due to the existence of the symmetry-breaking field. At criticality, 



the string structure of these states does not agree with the predictions of |42j - to solve 
this puzzle had been the original motivation for this work. Our numerical results suggest a 
scenario that may solve this apparent contradiction: the string type of the relevant excita- 
tions undergoes a number of changes as the field is switched on, and for large systems these 
reorganizations take place at particular values of the scaling parameter \i. Therefore, the 
string structure entering the massive scaling limit \i — » oo is that observed for large values 
of the elliptic nome q. 

Numerically, we have been able to identify the single particle states up to the fifth 
mass, and a number of two-particle states. Clearly, our results for the scaling functions are 
in complete agreement with both the analytic values of the mass ratios and with earlier 
numerical work on the Ising model in a magnetic field. As to the Bethe Ansatz solutions, 



apart from one exception, the string structure for large field is that conjectured in [42 . The 



one exception concerns the string type associated to the particle of mass m^. We found 
two states which contain this particle — one being the single particle state with conformal 
dimensions (3/2, 3/2), the other a two-particle state which also contains the lightest particle 
and has scaling dimensions (7/2,7/2) at criticality. In both cases, the observed string 



structures agree, but they differ from the one proposed in [42], compare Table IV. Of course, 



it cannot be ruled out completely that another change of string type appears at a larger value 
of the scaling parameter, but we found no trace of such a behavior for scaling parameters 
up to /i < 80. This clearly demands further clarification, and maybe the investigation of 
the dilute A 4 model can lead the way. In the scaling limit, this model is described by an E7 
theory of factorized scattering, but has so far eluded a TBA approach analogous to that of 
|42]| . In this case, we have been able to identify all seven single-particle states and observed 
interesting string solutions, details will be published soon p5] . 

A number of interesting questions arise in connection with the massive excitations of the 
critical dilute A 3 model. This is one property of the dilute A 3 model which distinguishes 
it from the proper Ising model in a magnetic field, since the critical Ising model does not 
show such excitations. Moreover, we suppose that this phenomenon is not particular to the 
specific model and will show up in the other dilute models as well. It would be interesting 
to understand the physical nature of these excitations, and to obtain analytic predictions 
for the observed mass ratios. Also, their dispersion relations have not been studied because 
we concentrated on momentum zero states throughout this work. 
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FIGURES 

FIG. 1. Effective adjacency diagram of the dilute A3 model. 
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FIG. 2. Sketch of a typical arrangement of the Bethe Ansatz roots Vj in the complex plane. 
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FIG. 3. The E 8 Dynkin diagram 
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FIG. 4. Low-energy part of conformal zero-momentum spectrum for the critical dilute A3 model. 
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FIG. 5. Solution of BAE for the second largest eigenvalue for varying nome q. Thick lines 
denote the real parts (normalized by — log(g)), thin lines the imaginary parts (in units of n) of the 
BA roots. 




-0.4 ■ 
-0.6- 
-0.8 ■ 
-1.0- 



0.0 ~ o!l 0.2 03 OA 0.5 0.6 oY o!<8 o'.9 " 1.0 

nome q 



20 



FIG. 6. Solution of BAE corresponding to the single-particle state with mass 771,4 m the scaling 
limit. The system size is N = 10. Thick lines denote the real parts (normalized by — log(g)), thin 
lines the imaginary parts (in units of tt) of the BA roots Vj. For each real part, the numbers in 
square brackets give the (approximate) imaginary parts (in units of tt/32) for large values of the 
scaling parameter /j. 
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FIG. 7. Same as Fig. |6], but for system size ./V = 20. For clarity, the numbers in square brackets 
have been omitted for all real roots. 




22 



FIG. 8. The phase function <p(v) ( |3.3[) for the solution of BAE corresponding to the sin- 
gle-particle state with mass 1714 in the scaling limit. The system size is N = 20, q = 1/5. Horizontal 
lines are drawn at integer values, crosses at intersections with the graph of <p(v) denote the real 
roots, and vertical lines indicate the positions of the real parts of the other seven roots. 
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FIG. 9. Scaling function Fq(/j,) ( p. 12 ) for the largest eigenvalue obtained from systems of size 
N = 50, 75, and 100. Here, the logarithm of F (ji) is plotted against // /15 = q 8 / 15 N ( pTTo| ) 
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FIG. 10. Scaling functions Fj(/j,) ( |2.11| ) for the excitations obtained from systems of size 
N = 100. 
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FIG. 11. Ratios of scaling functions Rj{jj) ( |2,13 ). The individual data points shown stem from 
systems of size N = 50, 75, and 100. Subscripts a, b, c distinguish different solutions with identical 
critical scaling dimensions. The Es mass ratios are indicated. 




26 



TABLES 



TABLE I. The nine thermodynamically significant string types of Ref. [42] and the correspond- 



mg masses. 


mass 


t 




£ (t) 


A(*)/5 


u 


n 
u 


1 
1 


n 
u 




mi 


1 


2 


1 


(-1,1) 


m 2 


7 


4 





(-3,-1,1,3) 


m 3 


2 


4 





(-4,-2,2,4) 


1JI4 


8 


5 


1 


(-12,-8,0,8,12) 


m 5 


3 


6 


1 


(-7,-5,-1,1,5,7) 


m 6 


6 


7 


1 


(-14,-6,-2,0,2,6,14) 


m-j 


4 


8 





(-10,-8,-4,-2,2,4,8,10) 


m 8 


5 


10 


1 


(-13,-11,-7,-5,-1,1,5,7,11,13) 



TABLE II. Finite-size approximants for central charge c and smallest scaling dimensions Xj. 



N 


c 




Z2 


%3 


X4 


X5 




x 7 


X8 


X9a 


x 9b,c 




10 


0.499 681 


0.125 080 


1.014 498 


2.196 


3.148 


4.272 


4.417 


5.449 


6.608 


6.771 


6.804 


7.942 


20 


0.499 920 


0.125 020 


1.003 543 


2.142 


3.034 


4.062 


4.191 


5.098 


6.142 


6.274 


6.274 


7.196 


50 


0.499 987 


0.125 003 


1.000 563 


2.128 


3.005 


4.010 


4.135 


5.015 


6.022 


6.148 


6.148 


7.030 


75 


0.499 994 


0.125 001 


1.000 250 


2.126 


3.002 


4.004 


4.130 


5.007 


6.010 


6.135 


6.135 


7.013 


100 


0.499 997 


0.125 001 


1.000 141 


2.126 


3.001 


4.002 


4.128 


5.004 


6.005 


6.131 


6.131 


7.007 


00 


1/2 


1/8 


1 


2+1/8 


3 


4 


4+1/8 


5 


6 


6+1/8 


6+1/8 


7 



TABLE III. Example of Bethe Ansatz solutions at criticality in sector i = 1 which differ from 
the ground-state solution for system size iV = 4 by a number of "massive" two-strings. Here, 
a = in/32. 



N = 4 


N = 6 


N = 8 


N = 10 


N = 12 


N = 14 


1.977 471 499 


1.980 014 437 


1.982 321219 


1.984 408 860 


1.986312 526 


1.988 062 914 


0.575 364 766 


0.577 361036 


0.579176 248 


0.580 821977 


0.582 324 931 


0.583 708 690 


-0.105 369 741 


-0.105 850 465 


-0.106 289 271 


-0.106 688 316 


-0.107 053 693 


-0.107 390 875 


-0.886 306 204 


-0.888 679 105 


-0.890 833 363 


-0.892 784 131 


-0.894 563 861 


-0.896 200 997 




-0.000 000 072 ± 
10.999 999 764 a 


0.020 936 097 + 
10.999 999 615cr 


0.035 854 325 ± 
10.999 999 489 a 


0.047 525 711 ± 
10.999 999 373(7 


0.057153 974 + 
10.999 999 263 a 






-0.020 936 294 ± 
10.999 999 730 a 


-0.000 000132 ± 
10.999 999 570(7 

-0.035 854 556 ± 
10.999 999 719cr 


0.014 769 565 ± 
10.999 999 459(7 

-0.014 769 869 + 
10.999 999 564(7 

-0.047 525 969 + 
10.999 999 712cr 


0.026 242 969 ± 
10.999 999 320(7 

0.000 000171 ± 
10.999 999 444ct 

-0.026 243 304 + 
10.999 999 563 cr 

-0.057154 256 + 
10.999 999 706(7 
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TABLE IV. Numerically observed string solutions for systems of size N = 100 for nome q = 
(fj, = 0) and q = 7/1000 (p ~ 39.4). Only the values of the non-real roots Vj are given. Here, 
g = — log(g) = — iiTT and a = iir/32. For those roots whose real part is very close to g/2 or whose 
imaginary part is almost an integer times <r, such that the number of digits given is not sufficient 
to see the slight differences, the superscripts + and — indicate on which side the actual data lie. 
We also include ranges of q values within which changes of the string patterns are observed. The 
three solutions which correspond to a critical scaling dimension of 6 + 1/8 are labeled by subscripts 



a, b, and c, where the latter two yield identical eigenvalues for finite systems. 



(A 


+ r, A + 


r) 


mass 


non-real roots at g = 


non-real roots at q = 7/1000 


pattern change 




(0,0) 

















(— — ) 

\ 16 ' 16' 




mi 


5.277 011 + 16 cr 


0.500 000+0 + 11.000 000"C7 


0.000 138 
0.000 161 


-0.000139 
-0.000170 




^2 ' 2> 




m,2 


-4.123 856+ 7.024170 c- 
— i.uoz (yu -p lu cr 


0.500 000" £>+ 4.999 988 cr 
n ^nn nnn— n 4- 1 ^ nnn ni 9 n 

U.OUUUUU ZH 1U.UUUU1Z (7 


0.000 560 


-0.000 640 






hi) 


m 3 


-3.791098 + 16(7 
3.719 808 + 11.045 724 ct 


0.500 000+ Q± 12.011549 a 
0.500000+0 + 9.988451 cr 


n nnn zL4n 

U. UUU 

0.000 580 


— n nnn 4^n 

U. UUU 1UU 

-0.000 640 




+ 1,| + 1) 


iri4 


-3.517 907 + 14.671118(7 
-3.458 496+ 5.430 532(7 
3.981 241 + 16 cr 


0.485 956 q ± 0.992132 cr 
0.486 904 £ ± 7.000 002 a 
0.486 904 g ± 12.999 998 cr 
0.500 000+0 + 16(7 


0.000 340 
0.002 900 


-0.000 350 
-0.002 950 




(2,2) 




2 mi 


-3.245 499 ± 10.982 954 a 
3.237651 ± 10.937395 a 


0.482 282 q ± 11. 000 000^(7 
0.517718 ± 11. 000000^(7 








+ 2 >Te~ 


1-2) 


mi + m2 


-3.282 380 + 10.815 346 a 
3.304 607 + 10.953 718 a 
4.238 695 + 16cr 


0.500 000+(?± 11.000 000-cr 
0.500 002 e ± 15.000 054 cr 
0.500 002 q ± 4.999 946 a 


0.001 000 
0.001410 


-0.001 100 
-0.001420 




+ 2,1 + 

(3,3) 


2) 


ms 
2 mi 


-3.356 336 ± 12.564 943 a 
-3.326 604+ 8.981386 (7 
3.297 948 + 16(7 

-2.875 085 ± 11.000 662 a 
2.875191 ± 11.000 554 a 


0.500 000^0+ 8.998 677cr 
0.500 000^0+ 11.000 000-C7 
0.500 000- ± 13.001 324 a 

0.455 576 ± 11.000 000+ a 
0.544424 ± 11. 000000+ (7 


0.000 800 
0.001275 


-0.000 900 
-0.001475 




+ 3,^- 


h3) a 


mi + m 3 


-2.890111 + 11.001019 (7 
2.889 954 ± 11.000 849 a 
4.741 946 + 16 cr 


0.500 000-0 + 11.000 000+ ct 
0.500 000+0+ 11.928013 cr 
0.500 000+0+ 10.071 987cr 


0.000 430 
0.000 520 
0.003 750 


-0.000 440 
-0.000 550 
-0.004100 






h3) 6 


mi + m2 


-3.555 509+ 8.138 452 (7 
-3.317 682 + 16(7 

2.889 954 + 11.000 849 cr 

4.741946 + 16cr 


0.471689 ± 11.000 000+ ct 
0.517 536 ± 5.000 047 a 
0.517 536 ± 14.999 953 a 


0.000 370 


-0.000 380 




+ 3,^- 


h3) c 


mi + m,2 


-2.890111 + 11.001019 (7 
3.091673+ 5.267 222 cr 
3.097 274 + 14.605 702 cr 


0.482 461 ± 5.000 047 a 
0.482 461 ± 14.999 953 cr 
0.528311 ± 11. 000000+ cr 






(1 


+ 3, |+3) 


mi + ni4 


-3.447863+ 6.844 587(7 
-3.399 457+ 16cr 
-3.048 514 ± 11.049 729 a 
3.054 397+ 11.080 356 cr 


0.485 373 ± 0.964 059 cr 
0.486 904 ± 7.000 017 cr 
0.486 904 ± 12.999 983 cr 
0.500 000+0 + 11. 000 000+cr 


0.003 800 


-0.003 850 



0.500 000+0 + 16(7 
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TABLE V. Ratios of scaled gaps for three values of the scaling parameter fx. 





N 


Ri 


R 2 


Rs 


Ra 


Rs 


Re 


R7 


RSa 


R»b,c 


R9 


10 


10 
20 
50 
75 
100 


1.618 749 
1.618 037 
1.618 024 
1.618 024 
1.618 025 


1.983 943 
1.983 923 
1.983 977 
1.983 984 
1.983 990 


2.403 558 
2.403 649 
2.403 797 
2.403 815 
2.403 834 


2.083 461 
2.081017 
2.079 804 
2.079 550 
2.079 551 


2.616910 
2.617335 
2.617468 
2.617482 
2.617 497 


2.910107 
2.930 942 
2.934 096 
2.934 583 
2.934 617 


2.317908 
2.320 478 
2.314 580 
2.313 631 
2.313 523 


3.123 561 
3.086 893 
3.078 762 
3.077 950 
3.077 537 


2.811391 
2.781193 
2.773 762 
2.772 773 
2.772 605 


3.360 590 
3.369 241 

3.375 544 

3.376 215 
3.376 509 


20 


10 
20 
50 
75 
100 


1.623 692 
1.618144 
1.618 034 
1.618 034 
1.618 034 


1.989 446 
1.987688 
1.987 597 
1.987591 
1.987 589 


2.412 441 
2.405 009 
2.404 866 
2.404 865 
2.404 865 


2.017 536 

2.025 316 

2.026 212 
2.026 286 
2.026 315 


2.623 691 
2.618144 
2.618 034 
2.618 034 
2.618 034 


2.958 322 
2.954 218 
2.954 289 
2.954 306 
2.954 311 


2.085 027 
2.129 306 
2.132 554 
2.132 780 
2.132 871 


3.006 016 
3.006 016 
3.003 342 
3.003 242 
3.003 212 


2.671041 
2.669 638 
2.669 034 
2.668 960 
2.668 941 


3.411712 
3.404 917 
3.404 805 
3.404 808 
3.404 808 


30 


10 
20 
50 
75 
100 


1.634 823 
1.618 424 
1.618 036 
1.618 034 
1.618 034 


1.992 979 
1.988 671 
1.988 497 
1.988 489 
1.988 486 


2.428 402 
2.405 376 
2.404 869 
2.404 867 
2.404 867 


2.003 678 
2.012146 
2.013 483 
2.013 604 
2.013 644 


2.634 823 
2.618 424 
2.618 036 
2.618 034 
2.618 034 


2.970 263 
2.956 329 
2.956 009 
2.956 007 
2.956 007 


2.021663 
2.071 838 
2.078 248 
2.078 742 
2.078 950 


2.994 370 
2.992 915 
2.992 952 
2.992 957 
2.992 958 


2.649 653 

2.644 399 
2.645130 

2.645 220 
2.645 248 


3.428 399 
3.405 376 
3.404 869 
3.404 867 
3.404 867 


oc 


oo 


1.618 034 


1.989 044 


2.404 867 


2.000 000 


2.618 034 


2.956 295 


2.000 000 


2.989 044 


2.618 034 


3.404 867 
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